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We derive the multiparticle-correlation expansion of the excess entropy 
of classical particles in the canonical ensemble using a new approach that 
elucidates the rationale behind each term in the expansion. This formula 
provides the theoretical framework for an entropy-based ordering criterion 
that has been already tested for a variety of model fluids and thermody- 
namic phenomena. In view of future investigations of the phase diagram 
of colloidal mixtures, we detail in this paper the case of a two-component 
system of spherical and rod-like particles and discuss the symmetries un- 
derlying both the self and distinct pair-distribution functions under various 
geometrical constraints. 
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1 Introduction 



One-phase criteria have been often introduced to estimate the location of the 
phase boundaries of the liquid and solid phases [1]. These empirical rules are 
usually rather specific. Nevertheless, they can be quite useful when it is not easy 
or straightforward to evaluate the free energies of the competing phases. Among 
such rules, one formerly proposed by Giaquinta and coworkers [2] has proved 
to be a rather general and flexible tool that can be successfully applied to a 
variety of models under different structural and thermodynamic conditions. This 
rule actually qualifies as an ordering criterion that can be readily implemented 
on the basis of the properties of the more disordered phase. The theoretical 
framework is provided by the multiparticle-correlation expansion (MPCE) of the 
excess entropy of a classical fluid that was originally derived by H. S. Green in 
the canonical ensemble [3], and was later extended by Nettleton and M. S. Green 
to an open system [4]. According to this formula, the entropy can be written 
as an infinite series whose n-th term is the contribution associated with density 
correlations involving n-particle multiplets. In short, the criterion states that the 
overall contribution to the entropy of a fluid arising from the spatial correlations 
involving more than two particles - a quantity called "residual multi-particle 
entropy" (RMPE) - changes sign concurrently with the local emergence of any 
new kind of structural organization in the system. The zero-RMPE criterion has 
been tested against such diverse thermodynamic phenomena as freezing [2, 5, 6], 
fluid-fluid phase separation in hard-sphere mixtures [7], mesophase formation in 
model liquid crystals [8, 9], the Kosterlitz and Thouless metal-insulator transition 
in a two-dimensional Coulomb lattice gas [10], and, more recently, the density- 
maximum anomaly in liquid water [11]. 

In a previous paper, we gave an entirely new proof of the entropy MPCE in 
the canonical ensemble by exploiting a simple combinatorial identity [12]. The 
proof applies to both continuous and lattice systems. In this paper we reconsider 
this derivation from a different perspective which discloses the significance of each 
term in the expansion. Moreover, we outline an iterative method for building up 
the expansion term by term. This method is then extended to systems composed 
of two different species of particles. Indeed, the present physical motivation 
for developing this formalism, besides the self-standing interest in the formal 
derivation of an entropy MPCE for binary mixtures, is that of demonstrating 
the effectiveness and reliability of the zero-RMPE criterion also in the case of 
a colloidal mixture of spherical and rod-like particles. In particular, we plan 
to check whether the miscible, low-density phase of such a model undergoes a 
kind of structural instability - of the type indicated by the vanishing of the 
RMPE - towards a more ordered phase, be it lamellar (for low concentrations of 
the spheres) or immiscible (the bulk-phase-separated system), for values of the 
total packing fraction close to those independently ascertained by experiment 
and numerical simulation [13, 14]. Previous studies have already shown that 



1 



the RMPE of parallel spherocylinders vanishes at about the same density where 
the smectic order sets in [8]. We are confident that the effects on the phase 
behavior that are associated with the addition of a small quantity of spheres will 
be correctly accounted for by the RMPE. 

While deferring to a forthcoming paper the discussion of the numerical simu- 
lation study of a mixture of hard spheres and spherocylinders, we anticipate here 
the analysis of the symmetries owned by the three pair distribution functions 
that enter the calculation of the two-body entropy. 



2 One-component systems 

We start discussing the entropy MPCE of one- component systems in the canonical 
ensemble. This choice is by no means restrictive since, as Baranyai and Evans 
first pointed out, it is always possible to take advantage of the canonical sum 
rules of the correlation functions and then rearrange the entropy MPCE in an 
ensemble-invariant form [15]. 

Let H N = {Ri, . . . , Rat} and P N = {Pi, . . . , Pjy} be the set of all particle 
coordinates and momenta, respectively. The canonical partition function of the 
system can always be split into an ideal and an excess part, = Z^Zff, where: 

^4Q" and Zf? = ±Jd'R 1 ...d*R„e-™W. (2.1) 

In Eq. (2.1), (3 = l/(k B T), A = h/^2nmk B T, V is the volume, and V N {R N ) is 
an arbitrary potential-energy function (in the most general case, Vn is a sum of 
n-body terms with n — 1,2, ... ,N). The excess entropy = Sn — reads: 



S%_ r d 3 i?x . . . d 3 R N e -^( RiV ) e~^ Rjv ) 
k B = ' V N Z% n Z% 

Upon defining a set of N normalized distribution functions (DFs) as: 

, -/3V N (R N ) 

P( N \R N ) = 



ycx 
N 



d 3 R n+1 ...d 3 R N e -WR") 

7ex 

with the properties 



p( n )m nx f ^R n+1 ...d^R N („ 1 N iwo^ 

P ( R ) = J f/iv^ ^ ( n = l,...,N-l), (2.3) 



J d^___d^ p( „) (RB) = 1 J d^+l pCn+D^n+l) = pW (R «) s 

(2.4) 
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the ordinary n-body DFs can be expressed as: 



f {n \r n ) = ( E 6 3 (R il -v 1 )...6*(R in -v n ) 

Nil- ..in I 



(JV _ n) ,/(Ri-r 1 )...^(R«-r r 

N\ pW(r") 
" ( N -n)\ V- ' (2 ' 5) 

where the sum is carried out over all n-tuples of distinct particle labels. The 
reduced n-body DFs (for n = 2, . . . , N) read: 



f(n) (r 



71 \ 



/W(r 1 ).../W(r B ) 
I M fl " " ' ' ' P'" ) (r" 



- ftfl-^^Hr"), (2.6) 



a=l x 

and verify the property: 

/^^pW(R n+ 0^ +1 >(R" +1 ) = (i - <? (n) (R n ) , (2-7) 

which holds also for n = 1 if = 1. We note that P (1) = 1 and = N/V 
if no one-body term is present in Vn, i.e., if no external potential acts on the 
particles (for the ideal gas, P^ = 1 for all n). From now on, we shall adopt the 
shorthand notations Pi 2 ...„ = P^(R") and Pi 2 ... n — P^(R"). Moreover, any 
integral of the form V~ n J d 3 Ri . . . d 3 R n (• • •) is hereafter denoted as /(• • ■). 

Now, we show how to build up the MPCE of the entropy term by term. The 
strategy is to consider a progressively increasing number of particles in the system. 
For N = 1, the (adimensional) excess entropy is just S^ x /k B = — J Pi In Pi, which 
leads to a first-order approximation to the excess entropy of a A^-particle system, 
Sfj/kB ~ S$/k B = —N J PxlnPx (viz., each particle contributes to the entropy 
independently of the other particles). 

Next, we move to a system of two particles only, and write its excess entropy 
as s!p plus a remainder k B R 2 that is equal to: 



Ccx C 

rCo = 



(l) 



k 



B 



J P 12 In P 12 + 2 J P 1 In Px = - J P 12 In P 12 . (2.8) 



Equation (2.8) leads to a second-order approximation for S^f under the hypothesis 
that each distinct pair of particles contributes an equal two-body residual term 
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to the entropy: 

7 (2) 



-N J Pi In Pi - My P 12 \nP 12 . (2.9) 



Note that the approximation (2.9) is exact for N = 2, i.e., = iSf*. 

When there are three particles in the system, the excess entropy is the sum 
of and a remainder kgR?,: 

R 3 = = _y P 123 lnP 12 3 + ^) / P i2lnPi2, (2.10) 

which suggests a third-order approximation for in the form 

g = - JV /p 1 l„P 1 _P)/p 12ln P 12 -(g [/P 123 lnP 123 -g)/P 12 lnA 2 

(2.11) 

Again, note that £3^ = Sg*. Equation (2.11) reproduces the first three terms in 
the r.h.s. of Eq. (5.9) of Ref. [12], and one may suspect that the remaining terms 
will be similarly discovered by arguing for N — 4, 5, . . . as we did for N = 1, 2, 3. 
The related proof can be obtained by induction over N. In fact, we know our 
target [12]: 

cex 

°V , ,, ,.. ,. 

12...N 



= -N J P, In P, - J P12...N In Pi 

= -iv/p.lnP.-EHEt-ir-fj/^lnP,... (2.12) 



After taking 7 a = / P L .. a lnP L .. a and J a = / P L .. a lnP L .. a (for a = 1, . . . , N), we 
first show that Eq. (2.12) is correct or, equivalently, that: 

'^C)s ( - ir tK (2 - i3) 

To this aim, it is enough to observe that the coefficient of I a in the sum (2.13) is 
(cf. Eq. (5.10) of [12]): 

£M3CH?; ^ N (214) 

Equation (2.12) can be built up term by term using the same procedure as 
sketched above for N = 1,2,3. Upon defining for a given M (with 2 < M < 
N — 1): 

= -nji - e ( N ) E(-ir a f!V- (2- 15 ) 



k B n=2\ n J a=2 
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which is tantamount to truncating the sum (2.12) over n after the M-th term, 
we set S c m +1 = + k B R M+1 , where: 



M 



rm+i = -(jM + i-(M+i)j 1 ) + ^[ M+1 )±(-ir- a l n )i a 



„=2 \ 71 J a=2 \ a s 



<M+1 




= -E 1 (-l) M+1 - a f M+1 Va- (2.16) 
a=2 V fl / 

This result allows one to define a higher-order approximation to the entropy as: 



(M + l) q (M) / AT \ 



n=2 \ n / a=2 
M+l 



Os-.r-C;')'. 



-^i-e ; E(-ir a ; 4, (2.17) 



n=2 \ / a=2 \ 



which is exactly the same as (2.15), but for the quantity M + l which replaces 
M. 

The MPCE of the entropy remains formally the same if the particles possess 
further degrees of freedom besides those pertaining to the centre of mass. For 
instance, in the case of liquid-crystal molecules (i.e., elongated particles with 
cylindrical symmetry), there are two more degrees of freedom for each particle 
since two angular coordinates are needed to specify the orientation of the molecule 
in three-dimensional space. If the body z axis is taken to coincide with the 
molecular axis, we can represent the direction of a molecule by resorting to the 
Euler angles 9 and <f) (see Fig. 1). The third angle, ip, describes a rotation around 
the molecular axis and, as such, it is not relevant for the configuration of the 
molecule. Let £ = (R, 9, <p) be the five-dimensional vector of coordinates of an 
individual molecule. The interaction potential is then a function of £ N . 

One of the simplest reference models for a liquid crystal is a system of hard 
spherocylinders or, equivalently, rods that cannot approach each other beyond 
a given distance a (the spherocylinder diameter). Such particles show up-down 
symmetry. In this specific case, the potential Vn{C, n ) will also reflect this sym- 
metry, in that it must be invariant upon interchanging (9, 0) with (jr — 9, n + <fi). 

The rotational kinetic energy K rot of a massive rod, written in canonical 
coordinates, reads: 
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where / is the moment of inertia relative to any axis perpendicular to the rod 
and passing through its centre. It then easily follows that the ideal and excess 
partition functions are: 

where A = h/ v / 27t7FbT and d 5 £ = sin 9 d 3 R d9 d(f>. The factor sin 9 in the volume 
element originates from the Gaussian integral over P^. Because of this factor, the 
delta function of argument £ should be intended as follows: 

5 5 (£ - 6) = 5 3 (R - Rx) 5{d 7y 5(0 - 0i) . (2.20) 

sin v 

With this proviso, the formal definition of the DFs, given by Eq. (2.5), as well as 
the overall appearance of the entropy formula (see Eq. (2.12)), remain unchanged 
provided that we now interpret /(•••) as (AixV)~ n J d 5 £i . . . d 5 £„ (• • •) (the value 
of n is always implicit in the form of the integrand) . 

Let us now consider the properties of the reduced pair distribution function 
(PDF) in relation to the symmetries of the system. Assume that no external 
field is present and that the molecular interaction includes, besides the hard-core 
repulsion, at most a pair term whose strength depends on the distance between 
the centres of mass. In such a case, the PDF g^(^i,^2) will only depend on the 
relative position of the two molecules. In a reference frame Si where molecule 1 
is placed at the origin and oriented along the z axis, the position of molecule 2 
is thoroughly described by the three (spherical) coordinates of its centre of mass 
( r i2)^i2) and (pi 2 ) plus two more angles, 6> 12 and <f>i 2 , specifying the direction of 
its axis. However, the above description is redundant since the orientation of, 
say, the x axis of Si is still arbitrary and we can always arrange things in such 
a way that (p 12 = 0. Hence, will depend on four variables only, namely one 
distance (r 12 ) and three angles (#12, #12, and (f>i 2 ). 

An even simpler situation is that of an artificially constrained nematic fluid, 
namely a system of elongated particles whose axes are kept parallel to each other 
while their centres of mass are free to move. In this case, the angles #12 and 0i 2 
are no longer necessary, with the result that the PDF depends just on r i2 and 
$12 (or, equivalently, on p 12 = r 12 sin-i9 12 and z 12 = r i2 cos$i 2 ). 

Another interesting case is that of an inhomogeneous system of rod-like molecules 
confined by an impenetrable wall. This model is useful for investigating the wet- 
ting properties of a nematic fluid and the onset and growth of a smectic layer at 
the wall. If the strength of the wall-particle attraction depends just on the dis- 
tance z of the molecular centre of mass from the wall, the number density 
turns out to be a function of x — ( z , 0) (i-e., the axis direction is relevant even 
when the wall does not exert any attraction), while the PDF ^ 2 '(^,^) depends 
on xi and Xi- 
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3 Two-component systems 

In this section, we generalize the MPCE of the entropy to binary systems com- 
posed of two different kinds of particles. We shall proceed in two steps: we shall 
first use the iterative method outlined in the previous section as a guidance for 
conjecturing the complete formula from its first few terms. Then, we shall give a 
formal proof of this formula by induction over the total number of particles. 

The canonical partition function of a two-component system with N = Ni+N 2 
classical point particles generally reads Zn 1) n 2 — ^nI 1 ^n 2 2 ^Ni,n 2 j wnere the 
excess part has the form: 

= ±J d^Rd^Qe-M^" 1 '* 1 *) . (3.1) 

In Eq. (3.1), the potential-energy function is arbitrary. The excess entropy is 
given by the integral 

-f^ = ~j d»*fld»*«5— ^ ln^^ . (3.2) 

KB V J ^Ni,N 2 N\,N 2 



Upon defining 



-0V Nl , N2 (B. N i,Q N 2) 
N\,N 2 



p(n u n 2)mni Q n 2) f d^ 1+1 . . . d^ d 3 Q ra2+1 . . . d*Q„ a ; ^ 

{ ^ ) J yJVi-m Zf? uN2 

(with n 1 <N 1 ,n 2 <N 2 ,n 1 +n 2 = l,...,N-l), (3.3) 

we have the following properties: 

r d 3 i?i . . . d 3 R ni d 3 Qi . . . d 3 Q n2 p ( ni) „ 2 )/ Rni nn2 \ _ -, . 
J V n > V n * [ ^ ' ~ ' 

, j3d 

J x ni+l p(ni+l,n 2 )^pni+l Q™2^ _ p(ni, n 2 ) ^ pni Q«2"j. 

y d Qn 2 +l p(n 1 ,n 2 +l) Qn 2 +1^ _ p(ni,n 2 )^pn! Q n 2^J (3 4) 

When external fields are absent, P^ 1 ' ) = 1 and p( 0,1 ) = 1. For a binary ideal 
mixture (i.e., Vn u n 2 = 0), the P functions are all equal to 1. 
The DF of order (rii, n 2 ) is defined as: 

, , jVJ iVo' p{ni,n 2 )f ni n 2 \ 
f{ni,n 2 )( r n 1 n 2 \ Jt l- Jt 2- -» v 1 ? _j I to c\ 

1 { ' q ] (iVx-m)! (iV 2 -n 2 )! • ^ 
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While the definition of the self reduced DFs g (ni '^ (for m > 1) and g {0 ^ (for 
n 2 > 1) is strictly analogous to that given for a one-component system, the 
distinct reduced DFs are defined, for ri\,n 2 > 1, as: 

5 (ni '' l2) (r ni ,q n2 ) = II (i - -j^-J II (l - P (ni ' n2) (r ni ,q n2 ) , (3.6) 

where 

p(ni,n2)( r m r.^ 1 ) 

p(ni,n 2 )( ri! n 2 \ — V ' " ' (O J) 

I »<1 J -p(i I o) (ri) ...p(i,o) (rni) p(o,i) (qi) ...p(o,i) (qn2) - ^'J 
For ni + n 2 > 1, a property analogous to that expressed in Eq. (2.7) holds: 



/ %^P (1 '°Vn 1+ i)<? (ni+1 '" 2 V 1+1 ,q" 2 ) = (l - ^) S ( 

(3.8) 

plus a similar identity involving g( ra i> ra 2+i) . 

/,From now on, a short notation is adopted where P ab stands for p( a > 6 )( r a , q fe ) 
and any integral of the form V~ m ~ n2 J d 3ni r d 3n2 g (• • •) is simply denoted as 
/(•••)■ Moreover, for future convenience, we take I a b = / P a blnP afe and J ab = 
/ P a felnP a b (with a < Ni, b < N 2 , a + b > 1). For instance, the excess entropy 
(3.2) can be rewritten as 



Ccx 

^1,^2 . ; p i p 

Pvijv 2 inrjVjjVj 



= — iVi | P 10 In Pio — N 2 J Poi lnPoi - | Pjv lJV2 lnP^ , (3.9) 

or JiV!Ar 2 = NiJi + N 2 J m + In 1 n 2 - 

We now move on to determine a MPCE for the entropy of a mixture. To 
this aim, we consider increasing particle numbers in the system, starting from 
N 1 + N 2 = 1 (i.e., only one particle, either of type 1 or 2, is present). By 
reasoning in the usual way, we immediately obtain a first-order approximation to 
the excess entropy in the form 



^N U N 2 

k B 



-iVx J P 10 In Pio - N 2 J Pqi In P 01 . (3.10) 



Then, we analyze the three cases with Ni + N 2 = 2. For N± — 2 and N 2 = 
(or the other way round), things are the same as for a one-component system, 
i.e., S^q = S^l + k B R 2 ,o, with P 2 ,o = -/ P2ohiP 20 . Instead, for Ni = N 2 = 1: 

R hl = ^ ^ = -y P n lnPn + y Pi lnPio + J P ilnP i = - J Pu InPn, 

(3.H) 
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thus leading to the following second-order approximation for S^f Ni : 



S 



(2) 

N U N 2 



k, 



-N x J P 10 In P 10 -N 2 J Poi In P i 



J P 20 In P 20 - N,N 2 J P n In P n - J 



Pq2 In Pq2 ■ 



(3.12) 

(3) 

For Ni + N 2 = 3, we just reproduce below the expression of S^ N2 , which 
results from carefully considering the implied four cases: 



o(3) 

k B 

JVi 
2 

3 



= -JVi J P w \nP w -N 2 J PoilnPoi 
P 2 o 



2 



P 30 In P 30 

y P12 lnPi2 - 2^ P n lnP n - ^ 



y P 2 olnP 20 - ( 
lnPn 



N 2 
2 

2 



P02 In P 02 



P02 In P02 



y p 20 inp 20 -2 y. 



P21 lnP 2 i - / P 20 lnP 20 -2 P n lnP n 



P 03 In P 03 - 



P02 In P02 
(3.13) 



On the basis of the above structure, we conjecture the following general for- 
mula: 

OCX 

Ni,N 2 



Ni J P10 In P w -N 2 J Poi In P i 

'NA (N 2 
v ni / U2. 



E 

ni + n 2 > 2 
ni < JVi , n 2 < ^2 



a + 6 > 2 
a < ni, b < n 2 



-a-b U A U 2 



Pah In P a b 



(3.14) 



which was preliminary checked against the form of N2 as independently got 
by our iterative method. Note that the double sum in Eq. (3.14) can also be 
arranged in such a way that all the terms with the same value of n = n\ + n 2 are 
gathered together: 

N1+N2 min{n,iVi} 



IN1N2 — E E 

n=2 n 1= max{n-Af2,0} 



'NA N 2 
.ni / \n - n 1/ 



E ( 

a + 6 > 2 
a < ni, 6 < n — n\ 



•Qn— a- 



n\ \ in — ni 



lab ■ 



(3.15) 

First of all, we prove that the Eq. (3.14) is an exact identity. Indeed, for given 
a and b (with a < N\, b < N 2 , a + b > 2), the coefficient of I a b on the r.h.s. of 
(3.14) is: 

'NA ( N 2 ^ 



E 

n-i + «2 > 2 
a < m < Ni , b < n 2 < N 2 



X -1 



■ ni+n 2 -a- 



-l "1 «2 



ni / V n 2/ 
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b lab - 



On account of Eq. (2.14), the number (3.16) is always zero but for the case (a, b) = 
(Ni,N 2 ), where it equals —1. Hence, Eq. (3.14) is correct. 

We next prove, by induction, that Eq. (3.14) also follows from our iterative 
method. We define for any M such that 2 < M < N — 1: 

q(M) 

k B 

M min{n,ATi} / Ar \ / » r \ / \ / 

- E E C „! 2 J E (-1)"-' 

n=2ni=max{n-W 2 ,0} \ V \ V a + b>2 

a < ni , b < n — n\ 

(3.17) 

Moreover, for each pair (M 1? M 2 ) satisfying Mi + M 2 = M+l, let us put S M : i M2 = 
^m\,m 2 + ^bRmi,m 2 i where: 

^Mi,m 2 = —(Jm 1 m 2 ~ MiJ w — M 2 J i) 

S v ' 

I M 1 M 2 

+ n? 2ni=m ax?-M 2 ,0}( n j( n - n i) a 2 ^ ( « ) ( 6 ) ^ 

a < ni , b < n — ni 

M+l min{n, Mi} 



-e e ; l: j e 

n=2 ni=max{n-M 2 ,0} \ ' 41 / \" ""V a + b>2 



a 



a < m, b < n — ni 

n=2 n 1= max{n— M 2 ,0} \ 1 / \ V a + b>2 , / \ 

a < ni, 6 < n — n\ 

Mi 



/MA/ M 2 \ f+i /^WM + l-m 

n 1= MTl-M 2 Ul J \ M + 1 " V a + b>2 ' ^ WV b , 

a < ni, 6 < M + 1 - ni 



'ab 



a + b > 2 
a < Mi, 6 < Af 2 



In the third line of the above equation, the identity (3.15) was used for Im 1 m 2 - In 
the end, the same structure of S^ N emerges for S^~^ (i.e., Eq. (3.17) where 
M + l replaces M), as long as we define, in close analogy with the one-component 
case: 

q (M+l) q (M) min{M+l,ATi} / at \ / AT \ 

^N U N 2 _ ^N U N 2 , ^ (Ml \ JV 2 \ , , 

fc S ^ Mi=max{M+l-JV 2 ,0} V M l/ V M + 1 " M l/ 
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This concludes our proof. 

Equations (2.12) and (3.14) express the canonical MPCE of the excess entropy 
in compact form for pure and mixed systems, respectively. In the literature, the 
expressions quoted for the first few terms of this expansion are far more involved 
because they are usually written in terms of the reduced DFs. We have already 
shown in [12] that the familiar form of the entropy formula for one-component 
systems in the canonical ensemble emerges when the P's are eliminated in favour 
of the (?'s through the reverse of Eq. (2.6). In doing so, a constant term can 
be extracted from the correlation integrals which, if absorbed into the ideal-gas 
part of the entropy, makes the latter equivalent to the entropy of the infinite- 
sized ideal-gas system (for the sake of clarity, this argument is reformulated in 
appendix A). 

Before providing a similar demonstration for mixed systems, let us see what 
happens to S^ N2 . Upon repeatedly using Eqs. (3.4) and (3.8), we obtain (with 



where all terms, including the two constants, are extensive, viz., each term scales 
in the thermodynamic limit linearly with either N or V. 

However, in a closed system the asymptotic value of the two-body seZ/ reduced 
DFs differs at large distances from 1 for 0(N~ l ) terms [16]. This makes the 
numerical evaluation of the integrals in (3.20) particularly harmful for a small 
system {i.e., sensitive to its boundary). To (partially) cure this problem, we can 
take advantage of the canonical sum rules for the DFs (see Eqs. (3.8)), adding 
(and subtracting) to each integral precisely the extensive term that makes the 
integrand of order N~ 2 at infinity: 



Pl = NJV and p 2 = N 2 /V): 




- = - Pl J d 3 rP 10 (r) lnP 10 (r) - p 2 J d 3 rP 01 (r) lnP 01 (r) 
ln iV^T"^ / d3rdVp io( r )^o(r')^o(r,r') In £ 20 (r, r') 




(3.20) 



- P1P2 
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1 N, 
— In ■ 



No 



No 



1, N 2 
— In ■ 



No-1 



(3.21) 



As a result, i) the contribution from the boundary now grows like V x V 2 ^ 3 N~ 2 oc 
iV~ 1//3 ; ii) the overall constant term outside of the integrals is of (9(1); and iii) the 
new integrals now conform to those in the grand-canonical-ensemble expansion. 
For general (Ni,N 2 ), the constant terms amount to 



x 




»"(:)(:) 



(3.22) 



For fixed a and b (with a + b > 2) the prefactor of each logarithm is the opposite 
of (3.16). Hence, the sum (3.22) equals l^AqlA^ 1 ) + ln^lA^ 2 ), leading to 
an alternative expression for the entropy MPCE: 



]NuN2 
kn 



= N, 



2 " MM?) 



+ N 2 



- - ln(p 2 A 3 ) 



E 

rti + 112 > 2 
«l < Nx , n 2 < N 2 




N ± J P 10 \nP 10 - N 2 Jp m hi F\ 



(3.23) 



We argue that, using the sum rules (3.8), it is always possible to "complete" 
the integrals / P a b In in such a way that the integrand behaves like N~ a ~ b at 
infinity. In a way analogous to the one-component case, this is accomplished by 
adding and subtracting a suitable constant (which is —Ni/2 — N 2 /2 to second 
order). If ensemble invariance must hold for the entropy MPCE, the constant 
needed is necessarily the opposite of the number that, added to 



3(3 X Mr— >(";)(> 

a < Til, 6 < Tl2 



max{0,a— 1} 

II " 

il=0 



AT,; 



max{0,6— 1} 

n 

«2=0 



1 - 



»2 
iV 2 



(3.24) 

gives back a (9(1) quantity as an outcome. 

If either ri\ or n 2 is zero, we are led to the one-component case, which is 
treated in appendix A. Otherwise, when both n\ and n 2 are non-zero, the sum 
in (3.24) yields 



(\ / \ f max {0,a— 1} 
:)(:){ s 

O < T»l, 6 < T12 



. max{0,6— 1} . 
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+ E(-!) 

a=0 



ni— a 



E (-!) ni 

a + 6 > 2 
«i + 1 < a < ni , b < ri2 

E (-!) ni 

a + 6 > 2 
a < ni , i'z + 1 < 6 < ri2 

Til —1 

E 

> 

U2 — 1 

I E 

12=1 



+» 2 -«4 «l »2 



+n 2 -a-6 "1 »2 



111 1 



ni 



a=«i+l 



ni-o n l 



In 1 - 



?2 



E (-i 

6=i 2 +l 



,n 2 -6 / n 2 



= 



(3.25) 



Hence, in order to produce an ensemble-invariant formula, a null integral cor- 
rection must be added to the canonical correlation integrals associated with the 
distinct DFs {n\ — n 2 — 1 was such a case, see Eqs. (3.20) and (3.21)). 
We reproduce hereafter the entropy MPCE in its final form: 



>Ni,N 2 



iVi 



No 



ln(p 2 A 3 ) 



No / Pni InPr 



E 

ni + ni > 2 
ni < Ni, n 2 < V 2 

- N x -N 2 . 



'NA (N 2 
nJ \n 2 . 



N, J Pi lnP 1( 

e 2 (->—(:)(:) 



a < ni, 6 < ri2 



(3.26) 



In the r.h.s. of the above formula, the first four terms are the only ones which 
survive in the absence of any interaction between particles and, as such, con- 
stitute the ideal contribution to the entropy. Instead, the last two terms arise 
from the resummation of the numbers (A. 7) for each species. When suitably ab- 
sorbed into the integrals / P a olng a0 and / P^hig^, these terms eventually make 
the canonical-ensemble entropy expansion look like the grand-canonical one. In 
particular, for a large (JVi, N 2 ^> 1) and homogeneous binary system, the entropy 
expansion starts as follows: 



S 



N lt N 2 



- In(Mi) 



No. 



5 



- ln(Mi) 



- \p\ J d 3 r dV [g 20 (r, r') ln# 20 (r, r') - g 20 (r, r') + 1] 

- P1P2 / d 3 r d V [g u (r, r') In g u (r, r') - g u (r, r') + 1] 

X -p\ J d 3 rdV[( ?0 2(r,r')ln( 702 (r,r / )-(?02(r,r / ) + l] + ...(3.27) 
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If the above expansion is truncated after its pair-correlation terms, this formula 
provides a low-density approximation for the entropy of the mixture. 

The MPCE of the entropy derived above applies to a binary mixture of inter- 
acting point particles. However, the formula does not change even if the particles 
possess further degrees of freedom. In fact, the same comments that we made at 
the end of section 2 still apply. A case of this sort is a mixture of spheres and sphe- 
rocylinders. We shall analyze in a future publication the phase diagram of this 
model system in terms of the RMPE. Here, we focus our attention on the general 
structure of the distinct two-body reduced DF, since the case of the sphere-sphere 
DF is obvious while the symmetries of the self DF for spherocylinders have been 
already discussed in section 2. 

In order to find the maximum number of independent scalar variables that 
intervene in the calculation of <7ii(£,q), it is convenient to work in a reference 
system where the sp hero cylinder (species 1) is centred at the origin and lies 
along the z axis. In this reference frame, the position of the sphere (species 2) 
can be parameterized in terms of two variables only, namely the length r 12 and 
colatitude $12 of the vector joining the two centres of mass (as usual, <^ 12 = for 
a convenient choice of the x axis). The same result is obviously obtained when 
viewing the situation from a reference system where the sphere is centred at the 
origin. In this case, the need for two further 6 12 and 0i 2 variables (see section 2) 
is only apparent, since the z and the x axis can be chosen in such a way that such 
two variables vanish altogether. If the axes of the spherocylinders are frozen and 
parallel to each other, no further simplification occurs, and the distinct two-body 
DF is again a function of r 12 and $12 . 

The numerical calculation of gn for a homogeneous mixture is carried out as 
follows. After taking £ = (r, 9, 0), x = q — r, and invoking homogeneity, we have: 

n 1 f I N\ N2 \ 

= ^7 (E/ dfld^e, - 9)5(^ - 0) J d 3 r5 3 (R t - r)5 3 (Q, - r - x)^ 

AT / N 2 \ 

1 'i;<S s (Q,--(Ri + x))) . (3.28) 



AttV 



\j=i 



Hence, p 2 <?n(£, q) is the average density of spheres in Ri+x when a spherocylin- 
der is centred in R 1 . An explicit formula is: 



AAA 2 (r 12 ,^ 12 ) 
p 2 ■ 2-KTi2 sin^i 2 Ari 2 A^i 2 



2n(r 12 , <? 12 ) - : ^JIT I! ^ > (3-29) 



AA/" 2 (r 12 , $ 12 ) being the number of spheres within a tiny spherical ring of volume 
27rr^ 2 sin^ 12 Ari 2 A'(9 12 , centred at the position specified by r 12 and i?i 2 . 
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4 Conclusions 



In this paper, we have outlined a constructive method for building up the entropy 
multiparticle-correlation expansion in the canonical ensemble, term by term, for 
both pure and mixed systems of classical particles. The aim of this effort is 
twofold: i) to unveil the hidden combinatorial structure behind the expansion; 
ii) to set the stage for an application of the entropy-based ordering criterion 
introduced by Giaquinta and coworkers to the phase diagram of a binary mixture 
of hard spheres and spherocylinders. In this respect, we have discussed here the 
general symmetries owned by the three pair distribution functions. A detailed 
analysis of the phase diagram of the model in the framework provided by the 
zero-RMPE criterion will be the object of a forthcoming publication. 
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A Ensemble invariance of the entropy MPCE 

In this appendix, an argument appearing in Ref. [12] is reproduced for the reader's 
convenience. This argument deals with the overall constant term that appears in 
Eq. (2.12) when we eliminate the normalized DFs in favour of the reduced DFs 
through Eq. (2.6). 

For one-component systems, this constant amounts to: 



E ( N ) E(- 1 ) n_a ( n ) ln {N ~ 1){N ~ 2) ' ' ' {N ~ a + 1] 



n=2 



n 



a=2 



N- 



a-1 



(A.l) 



For each a value, the prefactor of the respective logarithm is just the number 
(2.14). Hence, the sum (A.l) equals \n(N\N~ N ), thus yielding a new form of the 
entropy MPCE: 



k B 



\ ~ ln(M 3 ) 



N 



^ fihP '-S(»)S ( - ir 1")/ p ' Jn91 - 

(A.2) 

where p = N/V. 

We now show that any single term in the sum over n at (A.l) is extensive. In 
fact, we have first: 



n 



a=2 



k=i 



N 



31 



(A.3) 



We prove in appendix B that: 



E (-i) n " 

a=k+l 



n 



(-1) 



n—l—k 



n — 1 



Then, 



m^ y N - i)iN -^ iN - a+i} 



a=2 
jy\ n-l 

n 



E(-i 



k=l 



k 



N 



.n) I N 



n— 1 



E(-i: 



n—k 



n — 1 



n— 1 



(A.4) 



* + 5ivjE(-i)-* 



n — 1 



k=l 



k 2 + ... 



k=l 



(A.5) 



In appendix B, we also show that 

n-l 

fc=i V K 



0, for d = l,...,n-2 
— (n — 1)! , for d = n — 1 



(A.6) 
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In conclusion, we obtain: 



( N ) £(-p-°h , n £ - z3 - (» - " ± H „ __JL_ . (A . 7) 

W„2 V"/ 'V"" 1 n(n-l) 

As discussed in the main text, the integrals in Eq. (A. 2) cannot be easily com- 
puted numerically since, for finite N, the system boundary also contributes in a 
significant way. However, upon taking advantage of the canonical-ensemble sum 
rules for the reduced DFs (see Eqs. (2.7)), it should be always possible to make 
every integrand in Eq. (A. 2) sufficiently small at large distances. This is accom- 
plished by adding (and subtracting) a quantity equal to the number in (A. 7), with 
the result of leaving an overall 0(1) number outside of the integral. Furthermore, 
the new form of the integral can be made identical to the (so-called) fluctuation 
integral of the same order which appears in the grand-canonical-ensemble expan- 
sion. 

We show this explicitly for the third-order term in the entropy expansion, 
which, when expressed in terms of the reduced DFs, reads: 



Vl ln wkt^ - l p3 1 d3ri d3r2 dV3 p ^ ^ 111 - J ^~ 

3 J N(N - 2) 6 J 912913923 



In order to conform to the grand-canonical-ensemble expansion, we have to add 
the integral 

_I p 3 J d 3 ri d 3 r2 ^ ( _^ i23 + 3gugi3 _ + x) ( A _g) 

which, in view of the canonical-ensemble sum rules (2.4) and (2.7), is equal to 
—N/6. This number is exactly the same constant that must be subtracted to 

( \hi—, — - A.10 

V3/ N(N — 2) V ; 



in order to produce an 0(1) constant. As a caveat, we note that 



N 



-ip 3 J d 3 n d 3 r 2 d 3 r 3 P1P2P3 N123 + (1 - 2K)g 12 g 13 + (k - 2)g 12 + 1] = 

(A.H) 

for any real k, not simply —1. 
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B Two combinatorial identities 



In this appendix, the formulae (A. 4) and (A. 6) are proved by induction. 
First, we prove that, for any n > 2 and 1 < k < n — 1: 



E (-ir a 

a=k+l 



.a i 



n-l)-k 



'n-V 
k 



(B.l) 



Equation (B.l) is valid also for k = and any n > 1. We argue by induction over 
n. For n = 2 and fe = 1, the formula (B.l) is trivially correct. Then, assuming 
that the formula is correct for an arbitrary fixed n and all positive k < n, we 
calculate the l.h.s. of Eq. (B.l) for n + 1 and any 1 < k < n — 1 (the case k = n 
is obvious): 
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= (-1) 



n—k 



71 ' 



fc- l y 

(B.2) 



which is just the r.h.s. of Eq. (B.l), but for n + 1 which replaces n. 
Next, we show that, for n > 3: 
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— {n — 1)! , for d — n — 1 



while it is trivial to check that 
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Arguing inductively, let us suppose that Eq. (B.3) is valid for a given value n 
(checking this for n = 3 is immediate), and see what happens for n + 1: 
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- -&- ir i"*> j+ &- ir i"*ise) fcm+( - ir 

= -E(-ir*(": V+ e ff) fEt-W' 1 ; VI +(-d" 



/, I V k J m=0 \ m ) U=l V k 



(B.5) 



Focussing on this intermediate result, we distinguish three cases: 

i) For 1 < d < n — 1, the first term is zero as is the second for m > 0. As for 
the rest: 

E (-ir fe ( n ~ x ) + (-ir = - E(-i) (n_1)_ " ( n ~ x ) + 

= - E(-l) (n_1) "* (" 7 ^ = "(I - l)™" 1 = • (B.6) 
fc=0 V k J 

ii) For d = n — 1 instead, the first term is (n — 1)!, while the second is zero 
for < m < d. The rest equals to: 



n-l 



- 1)! + E(-ir fc ( n k l ) k n - 1 + E(-i)- fe ( n k l ) + (-1)" 



= E(-ir fc ( n /) + ( ~ 1)n = ' (R7) 

like in the previous case. 

iii) Finally, and given all the above results, for d = n we have: 

_g ( _ ir .(» - y + (J) g ( _ ir .(» - + g (») g ( _ ir , 

as we wanted to show. 
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FIGURE CAPTION 



1 : The Euler angles 6,<j), and ip: The axes x, y, and z form the laboratory 
reference frame, whereas x',y', and z' are parallel to the body set of axes. 
(, called nodal line, is the straight line perpendicular to z and z' . It can 
also be viewed as the intersection between the xy and x'y' planes. The 
z' axis is chosen so as to coincide with the symmetry axis of the molecule 
(represented in the picture as a rod). The Euler angles are in the ranges 
< 9 < 7T, < (f) < 2tt, and < ijj < 2n. 
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Figure 1: The Euler angles 6, 0, and ip: The axes x, y, and z form the laboratory 
reference frame, whereas x', y', and z' are parallel to the body set of axes. (, 
called nodal line, is the straight line perpendicular to z and z' . It can also be 
viewed as the intersection between the xy and x'y' planes. The z' axis is chosen 
so as to coincide with the symmetry axis of the molecule (represented in the 
picture as a rod). The Euler angles are in the ranges < # < 7r, < < 27r, and 
< ip < 2vr. 



23 



